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Abstract 

We review the properties of electron shuttles, i.e. nanoelectromechanical devices that transport 
electrons one-by-one by utilizing a combination of electronic and mechanical degrees of freedom. 
We focus on the extreme quantum limit, where the mechanical motion is quantized. We introduce 
the main theoretical tools needed for the analysis, e.g. generalized master equations and Wigner 
functions, and we outline the methods how the resulting large numerical problems can be handled. 
Illustrative results are given for current, noise, and full counting statistics for a number of model 
systems. Throughout the review we focus on the physics behind the various approximations, and 
some simple examples are given to illustrate the theoretical concepts. We also comment on the 
experimental situation. 

PACS numbers: 85.85.+j, 72.70.+m, 73.23.Hk, 73.63.-b 
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I. INTRODUCTION 



As advances in technology push the size of the electronic components towards the nanome- 
ter scale, the well-established technology of Micorelectromechanical Systems (MEMS) begins 
to acquire quantum features. This trend signals the birth of a new research field, Nanoelec- 
tromechanical Systems (NEMS). Despite of the infancy of NEMS, a large literature is already 
available on it, and for a broad overview the interested reader is referred to recent review 
article d 1 ^ 3 . The review at hand has a more restricted scope: it is devoted to a theoretical 
analysis of a very specific NEMS device, the electron shuttle. By presenting a detailed study 
of such an idealized model system we hope to be able to illustrate the basic physics and 
conceptual problems that need to be understood before a general theory of NEMS can be 
developed. 

The electron shuttle, originally introduced by Gorelik et al. in 1998 4 , consists of a mov- 
able nanoscopic grain which is coupled via tunnel barriers to source and drain electrodes. 
A subtle combination of quantum transport and the mechanical degrees of freedom play 
an essential role for the functionality of the device. The device can exhibit a "phase tran- 
sition": when a control parameter is tuned (this could be the damping of the oscillator), 
at a certain threshold value the system enters a new transport regime, the shuttle regime, 
where charge is transported in an orderly fashion (i.e. essentially a fixed number charges 
per mechanical oscillation cycle). (This phenomenon is particularly pronounced if the de- 
vice operates in the strong Coulomb blockade regime, when only one excess charge at a 
time is allowed in the movable part.) In the original suggestion of Ref. the motion of the 
movable grain was treated macroscopically, i.e., with Newton's equation of motion. In our 
group we have been asking questions of what happens when the movable body is so light 
that also its motion becomes quantized. For example, does the shuttle transition persist 
in the quantum regime? Another interesting question concerns the interpretation of mea- 
surements of NEMS devices, possibly exhibiting shuttling. As it turns out, a measurement 
of the stationary IV-characteristics does not always yield enough information to uniquely 
identify the underlying microscopic charge transport mechanism. A point in case is the Ceo 
SET experiment by Park et al.— where two alternative interpretations, namely incoherent 
phonon assisted tunneling^iM or shuttling,^*^ are plausible. The current noise provides an- 
other important characteristics, supplementary to the mean current . 11112113 The Fano factor 
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which characterizes the degree of correlation between charge transport events is a powerful 
diagnostic tool to distinguish between various transport mechanisms. Therefore, studies of 
the current noise in NEMS have become an active field of research . 14 i 15 i 16 i 17 i 18 i 19 i 20 Based on 
studies where the mechanical system is treated classically, one expects a giant enhancement 
of noise at the shuttling transition^, and we want to investigate whether this also occurs 
in the quantum case. One can go even further along this line: why stop at the current 
noise, which is the second cumulant of the full counting statics (FCS)? Surely FCS (which 
equals the probability P n (t) of n electrons being collected, say, in the right lead at time t) 
will contain even more detailed information. Consequently, FCS in mesoscopic devices has 
been a very active research field for some time now, and its importance and relevance have 
been underlined by a recent measurement of the third cumulant. Thus, FCS for NEMS is 
obviously an interesting issue. We are aware of a recent calculation for a classical, driven 
shuttle^, and at the end of this paper we illustrate some of our own very recent general- 
izations of these concepts to the quantum regime using a simple example; a full description 
will be given elsewhere 2 ^. 

We have investigated quantum shuttles in several recent papers , 19 i 22 i 23 i 24 i 25 i 26 where the 
full technical details can be found; the purpose of the present paper is to introduce some of 
the basic issues to a more general reader. Several other groups have also recently studied 



quantum shuttles, see e.g. Refs. 



27 



28.29 



3Cl The paper is organized as follows. In Section 



2 we introduce the models of quantum shuttles employed in our work. The total Hamilto- 
nian consisting of the "system" (both mechanical and electronic degrees of freedom of the 
quantum dot array), the leads, and a generic heat bath is used to illustrate the derivation 
of a Markovian generalized master equation (GME) which is the starting point of the the- 
oretical analysis. Along the way from the Hamiltonian to the generalized master equation 
we identify the necessary assumptions, and point out several issues of potential importance 
not addressed so far within the field of NEMS. 

In Section 3 we discuss the calculation of the current and the zero-frequency component 
of the current noise spectrum for a NEMS device described by a GME. We also give a 
qualitative discussion of the required numerical calculations, which are highly non-trivial due 
to the large dimensions of the involved matrices. Finally, some examples of our numerical 
results are displayed. Section 4 is devoted to an elementary discussion of the calculation of 
the full counting statistics, and a brief review of the experimental status. 
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II. THE MODEL 



At least two model systems have been considered in the literature. Gorelik et alA consid- 
ered a single movable quantum dot, while Armour and MacKinnon^ introduced a model of 
a three-dot array whose central dot is movable, see Fig^ Both of these systems are of in- 
trinsic interest and may find applications in real systems. The quantum system is assumed 
to be in the strong Coulomb blockade regime in which only none or one extra (spinless) 
electron in the whole system are allowed. Thus, in the one-dot case the electronic states are 
|0) and |1), while in the triple-dot case we have |0), \L), \C), and \R). The quantum system 
is coupled to two leads with a high bias applied between them. The bias is smaller than 
the charging energy for addition or removal of other electrons but otherwise it is the largest 
energy scale in the model. 

The moving dot interacts with its surroundings and its dissipative dynamics is described 
by the interaction with a generic heat bath. The Hamiltonian has the form 

H = H osc + H e \ + if e i_ie a ds(£) + -Pleads + Hb&th + -f^osc-bath (1) 



where 



H osc - — + -3- (la) 



describes the mechanical center-of-mass motion of the central dot as a one-dimensional 
harmonic oscillator with mass m and frequency Uq. We emphasize the importance of the 
nonlinear dependence of ff e i-ieads on the oscillator degree of freedom; the ^-dependence is 
often exponential which lies at the heart of the shuttling instability. We refer to the literature 
for explicit expressions for the other terms appearing in Eq.(£Q) 4,23,31 . The leads are held at 
different electrochemical potentials fiL,R whose difference gives the bias across the array (see 
Fig. (JTJ)). We assume that the tunneling densities of states are independent of energy. 
This is necessary for the first Markov approximation^ used later on, to hold. Further, as 
already implied above, we assume fii — > oo, hr — > — oo. These assumptions are necessary 
for the derivation of the Markovian dynamics of the movable system. Finally, the heat bath 
consisting of an infinite set of harmonic oscillators linearly coupled to the mechanical degree 
of freedom is described with the Caldeira-Leggett model^ 3 - in its Ohmic form. 
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A. Generalized Master Equation 



For the description of the model we use the language of quantum dissipative systems.— As 
the "system" (or "device" ) we take the electronic states plus the one-dimensional oscillator 
describing the center-of-mass motion of the central dot. The electronic leads and the heat 
bath interacting with the mechanical degree of freedom constitute the reservoirs. The task is 
to integrate out the degrees of freedom of the reservoirs to end up with an equation of motion 
for the system density operator. The derivation proceeds in two steps first integrating out 
the leads and then the heat bath in the weak coupling limit to get the desired GME for 
the system density operator. We remark that the assumed additivity of the two heat baths 
should be proven, however we are not aware of any such proof. 

As a specific example of an end result of this procedure, we give the GME governing the 
system density operator in the one-dot case, written here in a Liouvillean form: 

p(t) = Cp(t) = (£ coh + C driv + C daiap )p(t) (2) 

where the various superoperators C i: i= coh, driv, damp, are defined as 

£ C ohP = ^[-^osc + £qc\cq - eExclc , p] , (3) 
AmvP = — ^{cocle~~p - 2coe"pe"c + pe~~c cl) 

- ~2 ( c o c o e A P~ 2c o e A P e A c o + P e A c o c o ) , ( 4 ) 
C damp p = - t g[xA P ,p}]- : ^(N+l/2)[x,[x,p]) . (5) 

The physical meaning of certain of the terms is of particular interest. The term proportional 
to E is due to the electrostatic force on an occupied (i.e., charged) dot. The exponential 
dependence of the tunneling terms is clearly visible. The damping of the oscillator cannot be 
described entirely satisfactorily in the Markovian limit: one cannot simultaneously achieve 
translational invariance, positivity of the density matrix, and relaxation towards canoni- 
cal equilibrium. The form given above does not satisfy strict positivity; the quantitative 
consequences of this shortcoming turned out to be negligible^ 3 . To proceed further, one 
considers the electronic diagonal elements poo(t) = (0 |p(t)|0) and pn(t) = (1 |p(t)|l), where 
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|1) = c \0). These objects are still full density matrices in the phonon space and satisfy 
Poo(t) = -^[Hosc, Poo(t)] - -^-(e~~PooO) + poo(t)e~~) + T R e>p n (t)e* + C dSLmp p 00 (t) , 

1 xx T 2x 2x 

Pn(t) = ^[^osc - eEx,p n (t)} + T L e'~ poo(t)e~* - — (e~p n (t) + p n (t)e~) + C damp p n (t). 

(6) 

III. CURENT AND NOISE 

A. Analytical development 

The current through the system is given by 

/s tat = er L Tt omi (e-*f& t ) = eT R Tr osc (e^pff ) . (7) 

The trace is carried out over the oscillator basis and p^f = lim^oo p nn (t) . To find the 
stationary solution, one needs to truncate the oscillator basis at some suitably large value 
iV (in practice, N may reach 100 before convergence is achieved); since the resulting linear 
system has dimensions 2N 2 x 2N 2 (single-dot), or 10 iV 2 x 10 iV 2 (triple-dot), this becomes 
a critical issue. The method of solving this vast numerical task is commented on below. 

Let us next consider the calculation of noise. In quantum optics one often resorts to a 
result known as the Quantum Regression Theorem; this allows the calculation of any multi- 
time correlation function of system operators. Unfortunately, in NEMS this theorem can 
be applied only under very restricted circumstances: it can be used to calculate the noise 
(which is essentially the average of a product of current operators) within the three-dot 
device^ (because the current operators lie an d Irc operate within the "system"), while 
it is unapplicable for the single-dot case (and many other NEMS systems as well) because 
there the current operators involve operators belonging both to the baths (i.e., the electronic 
leads) and the system. A more general method is thus called for. 

In order to compute the noise spectrum, we follow the ideas of Gurvitz and Prager— , 
and introduce number- resolved density- matrices p\™\ where n is the number of electrons 
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tunnelled into the right lead by time t. Obviously, pu(t) = Y^nPu '(*)■ The p\™' obey 



» 



ih 



\H n 



\ e 

2 1 



PoO^ W] + ^damp Poo^ (^) 



Pn (t) = ^[Hosc - e£x,p^(t)] + £ damp pii ; (t) 
/// 



(8) 



with Pii (t) = 0. The mean current and the zero-frequency shot noise spectrum are given 
by^ 

d 



-f J. n ), 



n 

5(0) = ^! 



t— >oo 
2' 



(9) 
(10) 



where P n (i) = Tr osc [p^ (t) + Pi" (i)] are the probabilities of finding n electrons in the right 

lead by time t, i.e. precisely the objects needed for the FCS discussed below. We find 

— 

/ = Yl< n nP n {t) = r^Tr osc (e a p n (t)), i.e. one recovers the stationary current found above. 

— f i 

In a similar fashion, J2n n2p> riit) = TRTi osc [e a (2 np{l (t) + pii(t))], whose large-time 
asymptotics determines the shot noise according to (|10|). We have developed a generating 
function technique in Ref. Qto extract this large-time limit. Here we skip the technical de- 
tails; the upshot is that the zero-frequency noise, and thereby the Fano factor F = S(0)/2eI, 
can be expressed in terms of the pseudoinverse of the Liouvillean (and the static limit of the 
density matrix known already from the current calculation): 



2er 



R 



Tr 



e a 



QC~ X Q 



T R e^pff-e'' 




(11) 



J ii/ 



Here Q is a projection operator that projects away from the stationary state, for which the 
Liouvillean has the eigenvalue zero. The crucial point is that the pseudoinverse TZ of the 
Liouvillean, defined as QC^ 1 Q = 72. is tractable by similar numerical methods as used in 
the evaluation of the current. In our discussion of the full counting statistics given below, 
we analyze a toy-model to illustrate some properties of the pseudoinverse 71. 
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B. Comment on numerics; some results 



As mentioned above, the superoperator structure of the Liouville equation leads to large 
matrices of the order of N 2 x N 2 , where N is the number of low-energy states kept in the 
calculation. A further complication arises from the fact that the stationary limit corresponds 
to the zero eigenvalue of the Liouvillean, Cp = 0, forcing one to deal with singular matrices. 
The problems with the memory size can be circumvented by using iterative methods in 
which only LA for a given A is needed (N x iV numbers), and one avoids the storage of the 
full £ (N 2 x N 2 numbers). Using an iterative method raises the questions of convergence, 
and the speed of convergence. We found that the so-called Arnoldi iteration technique was 
sufficient for our purposes (see Appendix A in Ref. l25l ). provided that one uses a suitable 
preconditioning. While there exists a substantial empirical body of knowledge of how to 
carry out the preconditioning, we are not aware of a complete algorithm. As an example, 
when solving for the pseudoinverse TZ, we used the inverse of the "Sylvester part" Cq of L as 
preconditioned; finding the inverse of Co is a relatively fast procedure. When calculating 
the noise (which amounts to solving an equation of the type £x = b, where the vector b 
belongs in the range of C), we used the generalized minimum residual method (GMRes)^. 
Again, appropriate preconditioning was crucial. 

Figure El shows a set of current vs. damping curves for the single-dot shuttle. We 
draw attention to the following features. As damping is decreased, the current increases, 
approaching asymptotically the value / = 1/2tt ~ 0.16, i.e. one electron is transferred per 
cycle. In other words, this value of the current indicates that the shuttling transition has 
taken place, even in the quantum regime. The transition is not as sharp as found in the 
classical case-. For large values of damping the current is much smaller, and scales with the 
tunneling rate T. Very interestingly, we see a sharp increase in the current even without 
electric field, E — 0. Classically this does not happen, and therefore we interpret this 
cross-over of being due to quantum shot noise. 

We have in several occasions promoted the use of Wigner functions as an interpretative 
tool for the numerical results obtained for the stationary density matrix. The Wigner 
representation of the GME has also turned out to be a useful starting point for further 
analytic wor k 26 ' 29 ' 37 . These phase-space representations have a simple form in the classical 
limit: the Wigner representation of a regularly moving harmonic oscillator is an ellipse. On 
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the other hand, irregular motion under the influence of external noise gives rise to a Gaussian 
probability distribution centered at the origin. The charge resolved Wigner functions (n — 
corresponds to an empty dot, while n = 1 represents the occupied dot) are defined as 

w ^ p> - L <* - 1 "~ + 1 > exp (4>- < i2 > 

and some representative results are given in Fig. El As expected, as damping is decreased, 
the fuzzy central spot evolves into a ring. The finite thickness of the ring is due to thermal 
noise, the randomness of the charge transport processes, and the position-momentum un- 
certainty. The earmark of shuttling are the asymmetric, banana-shaped areas observed for 
weakest damping: here one observes a strong correlation between the occupancy of the dot, 
and the position and momentum. Thus, there is a large probability to have an occupied 
dot with a negative position coordinate and positive velocity (the dot has been filled in the 
neighborhood of the left contact), while there is a large probability of having an empty dot 
at positive x-values and negative velocity, i.e. on the return journey from the right contact. 
- As the system approaches the classical limit, the thickness of the Wigner ring shrinks; 
this is illustrated in Fig. As mentioned above, we expect that the noise of the quantum 
shuttle yields additional information about the nature of the charge transfer process. Figure 
E] shows some of our numerical results. In the top panel we see a relatively sharp cross-over 
if the tunneling length is larger than the length scale Xq = ^h~/(muj), i.e. as one aproaches 
the classical limit. In the bottom panel we show the Fano factor for the same parameters 
as in the top panel. One should note that a logarithmic scale is used, thus one observes (i) 
a huge enhancement at the tunneling cross-over, in particular for parameters approaching 
the semiclassical limit, and (ii) tiny Fano factors below the shuttling cross-over. Thus, even 
well in the quantum regime shuttling is a highly ordered charge transfer process. 

Figure El shows the phase-space plot for the parameter values indicated by the asterisk in 
Fig. El The interpretation is very suggestive: for these parameters tunneling and shuttling 
co- exist. Our numerics thus confirm the suggestion put forward in the analytical study of 
Ref. 29. The large value of the Fano factor can be understood as a slow switching process 
between the two possible current channels (tunneling and shuttling). In the next section we 
elaborate this point further. 
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IV. FULL COUNTING STATISTICS 



A. Calculational procedure 

We begin by introducing some notation. We recall that the Liouvillean, which is a non- 
hermitian operator, has a single eigenvalue equal to zero with p stat being the corresponding 
(normalized and unique) right eigenvector, we denote this eigenvector by |0)). The corre- 
sponding left eigenvector is the identity operator 1 which we denote by ((0|, and we have 
(( 1 )) = Tr(l''p stat ) = 1. The pair of eigenvectors allows us to define the complementary pro- 
jectors V = 1 0»«0 1 and Q=l-V obeying the relations VC = CP = and QCQ = C. We 
will also need the pseudoinverse of the Liouvillean 1Z = Q£~ 1 Q, which is well-defined, since 
the inversion is performed only in the subspace spanned by Q, where £ is regular. Rather 
than attempting to calculate the complete probability distribution P n (t) = Tr [/^(t)] di- 
rectly, it turns out to be easier to evaluate the cumulant generating function S(t, x), 



e s{t,x) 



J2 p n(t)e inx - (13) 



From S(t, x) we find the m'th cumulant of the charge distribution by taking the m-th 
derivative with respect to the counting field x at x = 0: 

am c 

«»-»(«) = s^U. (") 

and from the knowledge of all cumulants we can reconstruct P n (t). The cumulants of 
the current in the stationary limit t — > oo are given by the time derivative of the charge 
cumulants: 

«n> = 5«» m »(*)i^- (is) 

The first two current cumulants give the average current running through the system and 
the zero-frequency current noise, respectively. We have recently developed a systematic 
perturbation theory to calculate the higher cumulants^., and we just quote the final results: 



((I)) = «6|x|o», 

«/ 2 » = <<0|J|0>>-2((0|J^J|0>>, 

((/ 3 )) = «o|x|o» -3{(o\mj + jkx\q)) -6((o\m(mv -in)i\o)), (16) 
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where 



X 
J 



X' 



x- 
t 



(17) 
(18) 



The explicit expressions for the current superoperators X depend on the problem at hand; 
two examples are given in Ref. |22j (see also the analytic toy model discussed below). We 
emphasize the utility of these formulae: they imply that once the hard numerical problems 
with the stationary case are solved, the higher cumulants are essentially immediately ac- 



cessible. Using the systematic perturbation theory developed in Ref. |22j algebraic formulas 
for the fourth, fifth etc. cumulant can be generated with computer routines performing 
symbolic manipulations. In Fig. [7|we show numerical results for the first three cumlants for 
the single-dot shuttle. 



B. Toy model 

The general formulation for calculation of higher cumulants given above is very practical 
also in analytic calculations. Here we illustrate it by considering a toy model, for which the 
current and noise are very well-known, while the third cumulant cannot be considered as a 
standard result. Specifically, we consider a two level system, whose occupation probabilities 
obey the Master equation (p = (p Pi) T ) 

P = A> = \ VL Tr ]p (19) 

The relevance of this model to the shuttling dynamics is based on the fact that in the co- 
existence regime there are two "states" : the incoherent tunneling current channel (reflected 
by the central spot in the Wigner function description), and the orderly shuttling channel 
(reflected) by the ring, and that the system can be visualized as fluctuating between these 
two states. The bare two-level system, as discussed in this section, needs a slight modification 
to take into account that in the coexistence regime one deals with current channels and not 



states, and further discussion can be found in Ref. 



f. y 



22, see also Ref. 38 for a related model. 



The Liouvillean of Eq. (|T9*j) has the following right and left null vectors: 
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which satisfy £|0)) = = ((0|£ and ((0|0)) = 1, as is readily verified. We can immediately 
form the projectors V and Q: 

v = \oko\ = -^—[ TrTr ), q = i-v = -±-\ - ;-|. (2i) 

We also record the current superoperator: 




'-OR 




(22) 



Next, we need the eigenvector \u)) of L which does not belong to the null space (here, the 
situation is quite simple because the full problem is 2-dimensional with a 1-dimensional null 
space), and the associated eigenvalue \ v . We readily find 



K = -Fr + Tl)\ 14=1^; (^1 = fv+rl \ L Tr ) ' (23) 

and can thus compute the pseudoinverse TZ: 

n = QC-'Q = ±-\ V ))((P\ = - I f TL ~ Tr ) = yQ. (24) 

With these results at hand, the evaluation of various cumulants is reduced to simple matrix 
multiplications, and we quote the results 

«X» = =^£f- (25) 

J- L + t R 

« l2 » = (T*trlr m (26) 

= n - 2T%T L + 6ry| - 2TgT| + n 
(T R + v L y 



The first two results are well-known from the literature (see, e.g., Ref. |12J), however the 
third cumulant may be new. We believe that the above method of calculation can be very 
useful in extending several previous results, such as those obtained in by KieBlich et al.— . 



C. Experimental status 



To summarize the experimental status very shortly, we do not believe that the shuttling 
transition has been observed so far. Having stated this negative conclusion, we hasten to 

12 



point out that we think that the experimental situation is very promising, and that a crucial 
experiment could be just behind the corner. 

Two key experiments have been reported, which contain several important ingredients. 
In the experiment of Park et al. 5 a Cgo molecule was placed in a break junction, and the 
current- voltage characteristics showed clear indications of effects due to vibrational quanta. 
We believe that the observed values of the current are too small to be attributed to shuttling 
(recall that the shuttling current has the universal value of 1/2-71"), and that the system is in 
the tunneling limit. In the experiment of Erbe et a/.— the system was driven: indeed elec- 
trons were shuttled, but the experiment was not designed to observe the shuttling transition 
as a function of a control parameter. Our optimism is based on a number of new structures 
that are currently being explored in the literature, and in particular we find the structures of 
Scheible et al. 41 , where a soft silicon pillar forms the movable part, very promising. We also 
believe that the theoretical methods outlined above are suitable for modeling quantitatively 
the forthcoming experiments. 

V. CONCLUSION 

We have discussed at length some properties of a nanomechanical device, the quantum 
shuttle, which we believe can be an important component in future applications, for example 
measurements of very small displacements. While the quantitative results of the present 
work apply to very specific, and strongly idealized models, we believe that many of the 
phenomena we address are generic, and will be observed in near future. One of the central 
messages we want to pass is that the fluctuation properties of these devices contain a wealth 
of information, and that this information may be essential in identifying the key charge 
transfer processes in these devices. 
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► 

— Xq X Xq 

FIG. 1: Schematic picture of a three-dot system, introduced by Armour and MacKinnon 31 . The 
outer dots are fixed — the left one (L) at the position —xq and the right one (R) at xq, while the 
central one (C) can move (position x) in a harmonic confining potential. It also interacts with a 
heat bath causing damping and thermal noise. The outer dots whose respective energy levels are 
de-aligned by the device bias (eb) are coupled to the full or empty electronic reservoirs (leads), 
respectively. The current flows within the system due to tunneling between the left and central 
dot and the central and right dot. (Reproduced from Ref. I23L 
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FIG. 2: / — 7 curve. The damping dependence of the stationary current through the single-dot 
shuttle for different transfer rates and electric fields. Their values are d = 0.5:ro,r = 0.05/ia; 
(pluses; corresponds to Fig. |3J), d = 0.5xo,r = O.Olfaj (circles), d = 0.0, V = 0.05hto (asterisks), 
d = 0.0, r = 0.01k (crosses). Other parameters are A = xq,T = 0. The current is in units of eu> 



while 7 in huj. (Reproduced from Ref 



23) 
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FIG. 3: Phase space picture of the tunnelling-to-shuttling crossover. The respective rows show the 
Wigner distribution functions for the discharged (Woo); charged (Wn), and both (Wtot) states of 



the oscillator in the phase space (horizontal axis - coordinate in units of xq = y^h/muj, vertical axis 
- momentum in H/xq). The values of the parameters are: A = xo,T = 0, d = 0.5xo,T = 0.05hio. 
The values of 7 are in units of fko. The Wigner functions are normalized within each column. 
(Reproduced from Ref. l23h 
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FIG. 4: Transition from quantum limit to classical limit. The "quantum thickness" of the Wigner 
ring shrinks, and the radius increases, indicating larger maximal velocity, and oscillation amplitude. 
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FIG. 5: Current I and Fano factor F vs. damping 7. The 7-dependence of / (upper panel) and F 
on log scale (lower panel) for different transfer rates T and tunnelling lengths A. The parameters 
are A = xq,T = 0.01a; (full); A = xq,T = 0.05a; (long dashes); A = 2xq,F = 0.01a; (short dashes); 



A = 2xq,T = 0.05a; (dots) with xq = y h/muj. Other parameters are eE/ 



0.5x , T = 0. 



The current is in units of eu while 7 in units of uj. The asterisk defines the parameters of Wigner 
distributions in Fig. |H1 (Reproduced from Ref. 
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FIG. 6: Phase space picture of the shuttle around the transition where the shuttling and tunnelling 
regimes coexist. The respective rows show the Wigner distribution functions for the discharged 
(Woo), charged (Wn), and both (Wtot = Woo + Wu) states of the oscillator in the phase space 
(horizontal axis - coordinate in units of xq = y^h/muj, vertical axis - momentum in H/xq). The 
values of the parameters are: A = 2xq, eE/muJ 1 = 0.5xq, 7 = 0.029a;, T = 0.01a;, T = 0. 
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FIG. 7: The first three cumulants for the one-dot shuttle as a function of the damping 7. The 
parameters are A = 1.5xq, d = eE/muJ 2 = 0.5xo- 
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